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Abstract. Previous timing data for PSR B1620— 26 were consistent 
with a second companion mass ni2 anywhere in the range ~ 1CP 3 — 
IMq, i.e., from a Jupiter-type planet to a star. We present the latest 
timing parameters for the system, including a significant change in the 
projected semi-major axis of the inner binary, a marginal detection of 
the fourth time derivative of the pulse frequency, and the pulsar proper 
motion (which is in agreement with published values for the proper motion 
of M4), and use them to further constrain the mass m<i and the orbital 

parameters. Using the observed value of / , we obtain a one-parameter 
family of solutions, all with mi 10~ 2 Mq, i.e., excluding stellar masses. 

Varying / within its formal la error bar does not affect the mass range 
significantly. However, if we vary / within a 4<r error bar, we find that 
stellar-mass solutions are still possible. We also calculate the predicted 
rate of change of the projected semi-major axis of the inner binary and 
show that it agrees with the measured value. 



1. Introduction 



The millisecond radio pulsar PSR B1620— 26, in the globular cluster M4, has a 
low-mass binary companion (probably a white dwarf of mass m\ ~ 0.3 Mq for a 
pulsar mass m p = 1.35 Mq) in a 191 day low-eccentricity orbit. In addition, its 
unusually large frequency second and third derivatives indicate the presence of a 
second companion in a wider orbit around the inner binary (Backer et al. 1993; 
Thorsett et al. 1993; Michel 1994). Such a hierarchical triple configuration 
is expected to be produced quite easily in a dense globular cluster through 
dynamical interactions between binaries (Rasio et al. 1995; Sigurdsson 1995). 
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2. Latest Timing Data 

Timing data were taken at the VLA and Green Bank (140-ft), using standard 
hardware and techniques. Observations at the VLA were made at 1.6 GHz; ob- 
servations at Green Bank were made at a variety of frequencies between 390 MHz 
and 1.6 GHz. The data span is 1988 March to 1995 December. 

Timing parameters are reported in Tablet, and in most cases represent 
an incremental improvement over previously published values. We report for 
the first time a value for / , but emphasize that despite its high significance 
in this fit, covariance with unmodeled parameters (e.g., / ) may complicate 
its interpretation. An identical fit with a fifth order polynomial yields / = 
— 1.5±2.2x 10~ 40 s -5 and / = — l±6x 10 8 s -6 (the other frequency derivatives 
are unchanged). 

The globular cluster proper motion in RA and Dec is —9.7(7) and —12.4(7) 
mas/yr QCudworth fc Hansen 19"93| ), This is about 2a from our timing proper 
motion, a difference that we do not believe is significant. 

The optical proper motion has mag nitude \x = 2.42 x 10~ 15 rads -1 . This 
will produce a change in the projected semi-major axis x = aisini of x/x = 
/xsinj'/ tani, where i is the orbital inclination and j is the angle between the 
proper motion vector and the line of nodes (neither i nor j is known). If the 
observed x is due to proper motion, then sinj < 1 implies i < 15°, which in 
turn implies that the companion must have mass > 1.6Mq, i.e., it is a neu- 
tron star or black hole, not a white dwarf. The nearly circular orbit is then 
problematic, but could be explained if the companion were captured in a three 
or four body interaction which involved a collision, followed by circularization 
during a Thorne-Zytkow phase or by a disk of material from a disrupted star. 
In our opinion, however, it is more likely that the observed x has, instead, some 
contribution from another source, such as precession. 

At the most likely orbital inclination i = 60°, proper motion can account 
for only about 10% of the observed x. A more attractive explanation is that x 
is caused by precession of the binary orbit in the gravitational field of the same 
third body that has been invoked to explain the large spin-frequency derivatives. 
It can be shown that this "planetary" precession yields x/x = G1712F / '(Qir^) , 
where fii is the angular velocity of the inner orbit and r 2 is the distance to 
the second companion (assumed fixed). F is a complex function of unknown 
angles, but has a median value for random orientations of 0.12, a 50% range of 
0.033 - 0.33, and a 95% range of 0.0013 - 1.6. 

Assuming the inner binary has mass M = m p + mi = 1.7 M@ and a low- 
eccentricity outer orbit, and ascribing the entire x to precession, we can write 
F = (P2/3360 yr) 2 (l + M/m 2 ), where P2 is the period of the outer orbit, and 
calculate F for the three illustrative solutions of Rasio (1994). For mi = 80® and 
P 2 = 10 yr, F = 0.06, while for m 2 = O.8M and P 2 = 120 yr, F = 0.004, and 
for 777,2 = 1.4M Q and P 2 = 500 yr, F = 0.04. We conclude that the measurement 
of x is consistent with the presence of a third companion, and indeed is a strong 
indication that the triple hypothesis is correct. However, by itself it is only a poor 
discriminator between stellar and planetary solutions. With more information 
from spin frequency derivatives and from other orbital perturbations such as to 
and Pb, the measurement of x will help untangle the various unknown angles 
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that characterize the system. A complete analysis of these orbital perturbation 
effects will be presented elsewhere (Joshi & Rasio 1996). 

3. Orbital Parameters 

The best test to confirm the triple nature of this system would be a Keplerian 
fit to the timing data spanning more than one orbit of the outer body. However, 
since the inferred orbital period P2 is of order a century or more, we observe 
only a small portion of an orbit, effectively measuring successively higher or- 
der derivatives of the acceleration at a single point. If derivatives of the pulse 
frequency up to the fifth order were available, the system could in principle be 
solved completely (except for the usual inclination angle) for the orbital parame- 
ters and mass of the second companion. Here we use the values listed in Table 1 
for frequency derivatives up to the fourth order to obtain a one-parameter family 
of solutions. While computing the orbital parameters of the second companion, 
we treat the inner binary as a single object of mass M = 1.7 Mq. For fixed val- 
ues of the eccentricity e2 and the inclination 12, we solve the non- linear system 
of equations for the frequency derivatives (using the Newton-Raphson method) 
for the mass 777,2, the semi-major axis 02, the angle of periastron u>2 and the 
longitude from pericenter A2 (Joshi & Rasio 1996). 

Fig. 1 illustrates our "standard solution" , obtained using the observed value 

of /. We see that there are no solutions for 62 ^ 0.1. Hence a nearly circular 
orbit is ruled out. For 0.1 ^ e2 0.3 there are two solutions for each value 
of the eccentricity, and hence two possible values of 7772. In one solution, 7772 
approaches zero as e2 approaches ~ 0.3. However for 7772 SS 10~ 5 M@ the triple 
configuration becomes dynamically unstable. In addition, we can rule out these 
solutions because they have very short orbital periods (P2 SS 5yr) and this 
would have been detected already in the timing data. For the other solutions, 
7772 increases monotonically in the range ~ 10~ 3 — 10~ 2 M Q (Jupiter to brown- 
dwarf masses) as e2 increases from ~ 0.1 to 1. 

In Fig. 2 we show the results obtained on varying / within a 4<r error bar 
around the best fit value f m (given in Table 1). The results are shown for 
/ / J m = 1.0 (as in Fig. 1, solid line), 2.3 (long-dashed line), 0.001 (short- 
dashed line), —0.05 (dotted line) and —0.4 (dot-dashed line). We find that 

stellar-mass solutions for 7772 are still possible if —0.05 <J/ / / m <^ 0.001. A 
second companion of stellar mass would provide a natural explanation for the 
eccentricity of the inner binary in terms of secular perturbations (Rasio 1994) 
and would also be consistent with a preliminary identification of an optical 
counterpart for the system (Bailyn et al. 1994). 
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Figure 1. Our "standard solution," using the observed value of /. 
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Figure 2. The solutions for f = f m (solid line), for /= 2.3 f m (long- 
dashed line), for /= 0.001 f m (short-dashed line), for /= —0.05 f m 

(dotted line) and /= —0.4 f m (dot-dashed line). We assume here 
that the inclination angle %2 = 90°. 
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Right ascension (J2000.0) 
Declination (J2000.0) 
Proper motion RA (masyr -1 ) 
Proper motion Dec(masyr _1 ) 
Dispersion measure ( cm -3 pc) 

Spin period P (ms) 
Spin frequency / (Hz) 

/ (s- 3 ) 

7 (s- 4 ) 
7 r 5 ) 

Epoch of / (MJD) 



16 h 23 m 38!2228(6) 

-26°31'53 / /74(4) 

-16(3) 

-30(18) 

62.8627(8) 

11.075750892214(6) 
90.2873321847(5) 
-6.065(2) x 1(T 15 
1.897(4) x 10~ 23 
1.2(2) x 10~ 32 

-1.7(5) x 10- 40 
48365.0 



Projected semi-major axis x (s) 64.809478(8) 



Orbital period Pb (s) 
Eccentricity e 

Time of periastron T (MJD) 
Angle of periastron lo 
Mass function (M ) 

Advance of periastron oj (°yr _1 ) 
Pb 

e (s- 1 ) 
x 



16540653(6) 
0.0253151(3) 
48345.3771(3) 
117.1296(6) 
7.975 x 10~ 3 

(-2.9 ± 2.4) x 10~ 4 
(1.7 ±2.6) x 10" 9 
5(4) x 10~ 15 
6(1) x 10~ 13 



Table 1. Timing parameters of PSR B1620— 26. Position is relative 
to the JPL DE202 solar system ephemeris. Numbers in parentheses are 
uncertainties in the final digits quoted. NOTE: formal uncertainties 
are relative to model fit with above parameters. Covariances with 
unfit parameters (e.g., the fifth frequency derivative) may increase true 

uncertainties, particularly of / and / . 
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